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Abstract 



n 

Ph. In 111 we introduced an inflationary scenario, non-abelian gauge field inflation or gauge-flation 

for short, in which slow-roll inflation is driven by non-Abelian gauge field minimally coupled to 



gravity. We present a more detailed analysis, both numerical and analytical, of the gauge-flation. 
By studying the phase diagrams of the theory, we show that getting enough number of e-folds 



> 

£f~) ■ during a slow-roll inflation is fairly robust to the choice of initial gauge field values. In addition, 



we present a detailed analysis of the cosmic perturbation theory in gauge-flation which has many 
special and interesting features compared the standard scalar-driven inflationary models. The 
specific gauge-flation model we study in this paper has two parameters, a cutoff scale A and the 
gauge coupling g. Fitting our results with the current cosmological data fixes A ~ 10 H ~ 10 15 
GeV (H is the Hubble parameter) and g ~ 10~ 4 , which are in the natural range of parameters in 
generic particle physics beyond standard models. Our model also predicts a tensor-to-scalar ratio 
r > 0.05, in the range detectable by the Planck satellite. 
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I. INTRODUCTION 



The idea of inflationary cosmology was originally proposed to provide a possible resolution 
;o some of the theoretical problems of the bing-bang model for the early Universe cosmology 
2J1. However, with the advancement of the cosmoloeical observations and most notably 

n n 

the cosmic microwave background (CMB) observations [2j, |3[, the inflationary paradigm 
has received observational support and inflation is now considered an integral part of the 
standard model of cosmology with the following general picture. A patch of the early 
Universe which is a few Planck lengths in size under the gravitational effects of the matter 
present there undergoes a rapid (usually exponential) expansion, the inflationary period. 
The inflation ends while most of the energy content of the Universe is still concentrated in the 
field(s) driving inflation, the inflaton field(s). This energy should now be transferred to the 
other fields and particles, the (beyond) standard model particles, through the (p)reheating 
process. The rest of the picture is that of the standard hot big-bang scenario, with radiation 
dominated, matter dominated and finally the dark energy dominated era that we live in. 

In the absence of a direct observation for the primordial gravity waves, one of the main 
standing issues in inflation is what is the Hubble parameter during inflation H, or the 
energy density of the inflaton field(s). With the current observations, and within the slow- roll 
inflation scenario, the preferred scale is H < l(T 5 M pl , where M pl = (87iG N )-^ 2 = 2.43 x 10 18 
GeV is the reduced Planck mass. On the other hand, according to the lore in beyond 
standard particle physics models, the supersymmetric grand unified theories (SUSY GUTs) 
{4], the unification scale is around 10 16 GeV, suggesting that inflationary model building 
should be sought for within various corners of such models. If so, the SUSY GUT setting 
will also provide a natural arena for building the (p)reheating models. 

Almost all of inflationary models or model building ideas that appear in the literature use 
one or more scalar fields with a suitable potential to provide for the matter field inducing 
the inflationary expansion of the early Universe. The choice of scalar fields is made primar- 
ily because we work within the homogeneous and isotropic Friedmann- Robert son- Walker 
(FRW) cosmology and that turning on spinor or gauge fields generically violates these sym- 
metries. Another reason is that, from the model building viewpoint, turning on potential 
for the scalar fields is easier than for other fields, whose interactions are generically fixed 
by gauge symmetries or renormalizability conditions. Building inflationary models within 
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the SUSY GUTs then amounts to exploring various corners of the theory/model in search 
of fiat enough potential which supports successful slow-roll inflation, the flatness of which 
is respected by the loop and quantum corrections. Such models usually come under the 
D-term or F-term inflationary models 5]. 

Regardless of the details, non-Abelian gauge field theories are the widely accepted frame- 
work for building particle physics models, and, in particular, beyond standard models and 
GUTs. In view of the ubiquitous appearance of non-Abelian gauge fields, one may explore 
the idea of using gauge fields as inflaton fields, the fields which get nonzero background 
value during inflation and drive the inflationary dynamics. One of the main obstacles in this 
regard is the vector nature of the gauge fields and that turning them on in the background 
will spoil the rotation symmetry. 

A related scenario in which this problem was pointed out and addressed is "vector inflation 
6|." The idea in vector inflation, unlike ours, is to use vector fields and not gauge fields, 
as inflaton. In two possible ways were proposed to overcome the broken rotational 
invariance caused by the vector inflaton fields: (1) introduce a large number of vectors 
each assuming a random direction in the 3D space, such that on the average we recover 
a rotational invariant background; or, (2) introduce three orthogonal vector fields of the 
same value which act as the triad of the spatial part of the spacetime, the "triad method" 



[8j. The other important obstacle in the way of driving inflation by vector fields is the 
exponential, l/a(t) suppression of the massless vector fields in an inflationary background, 
causing inflation to end too fast. This problem has been overcome by adding nonminimal 

n n 

coupling to the gravity, usually a conformal mass type term [6|, |9[. To have a successful 
inflation, however, this is not enough and one should add quite nontrivial potentials for the 
vector field 0, Q, 0]. Dealing with vector fields, and not gauge fields, may bring instabilities 
in the theory: the longitudinal mode of the vector field which has a ghost type kinetic term 
and is not dynamical at tree level, in the absence of gauge invariance, can and will, become 
dynamical once quantum (loop) effects are taken into account. This latter will cause ghost 
instability, if we were studying the theory on a flat background. It has been argued that 



such instabilities can persist in the inflationary background too [10| ; see, however, ll| for a 
counter argument. In any case the instability issue of vector inflationary models seems not 
to be settled yet. 

In order not to face the above issue one should build a "gauge invariant vector inflation." 
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One can easily observe that it is not possible to get a successful inflation with some number 
of U(l) gauge fields. The other option is to consider non-Abelian gauge theories. The "triad 
method" mentioned above is naturally realized within the non-Abelian gauge symmetry 
setting, irrespective of the gauge group in question. We then face the second obstacle, the 
l/a(t) suppression. This may be achieved by ch ang ing the gravity theory, considering Yang- 
Mills action coupled to F(R) modified gravity 12], or considering Einstein gravity coupled 
to a generic (not necessarily Yang- Mills) gauge theory action. This latter is the idea we will 
explore in this work. We should stress that, as will become clear, our approach and that 



of 



12) are basically different. Using non-Abelian gauge fields has another advantage that, 



due to the presence of [A^, A u ] term in the gauge field strength it naturally leads to 
a "potential" term for the gauge fields which, upon a suitable choice of the gauge theory 
action, can be used to overcome the l/a(t) suppression problem mentioned above. 

In this work, we present a detailed discussion and analysis of gauge-flation, inflation 
driven by non-Abelian gauge fields, which we introduced in |l| . In section [Til we show 
how the rotation symmetry breaking can be compensated by the SU(2) (sub)group of the 
global part of non-Abelian gauge symmetry and how one can introduce a combination of 
the gauge field components which effectively behaves as spacetime scalar field (on the FRW 
background); and that there is a consistent truncation from the classical phase space of the 
non-Abelian gauge theory to the sector which only involves this scalar field. 

Setting the stage, in section IHIl we choose a specific action for the gauge theory that 
is Yang-Mills plus a specific F A term which can come from specific (one) loop corrections 
to the gauge theory. In this work, however, we adopt a phenomenological viewpoint and 
choose this specific F 4 term primarily for the purpose of inflationary model building. The 
important point of providing field theoretical justifications for this F A term will be briefly 
discussed in the discussion section and dealt with in more detail in an upcoming publication. 
Our model has hence two parameters, the gauge coupling g and the ceofficient of this specific 
F A term k. These two parameters will be determined only by focusing on the considerations 
coming from cosmological observations. In this section we present an analytic study of 
the inflationary dynamics of our gauge-flation model and show that the model allows for a 
successful slow-roll inflationary period which leads to enough number of e-folds. In section 
IIVI we present the diagrams and graphs for the numerical analysis of the gauge-flation 
model. Our numerical analysis reveals that the classical slow-roll inflationary trajectory is 
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fairly robust to the choice of initial conditions. 

Having studied the classical inflationary dynamics, in section [V] we turn to the question 
of cosmic perturbation theory in the gauge-nation. Because of the existence of other compo- 
nents of the gauge fields which has been turned off in the classical inflationary background, 
the situation here is considerably different than the standard cosmic perturbation theory 
developed in the literature. We hence first develop the cosmic perturbation theory for our 
model, discuss its subtleties and novelties; we discuss the scalar, vector and tensor perturba- 
tions, their power spectra and the spectral tilts. In section IVI} after completing the analysis 
of the model, we confront our model with the available cosmological and CMB data. We 
show that indeed it is possible to get a successful inflationary model with the gauge-flation 
setup. In the last section we summarize our results and make concluding remarks. In three 
appendices we have gathered some technical details of the cosmic perturbation theory. 

II. THE SETUP 

In this section we first demonstrate how the rotation symmetry is retained in the gauge- 
flation and then discuss truncation to the scalar sector. Here we will consider an SU(2) 
gauge theory with gauge fields A ^ where a = 1,2,3 label the gauge algebra indices and 
/i = 0,1,2,3 the spacetime indices, the temporal components will be denoted by Aq and 
the spatial components by A\. Although we focus on the SU(2) gauge theory, our analysis 
holds for any non-Abelian gauge group G, as any non-Abelian group always has an SU (2) 
subgroup. 

We will consider gauge- and Lorentz-invariant theories, where the gravity part is the usual 
Einstein-Hilbert action and the Lagrangian of the gauge theory part, which is minimally 
coupled to gravity, is of the form £ = £j(F^ v ; g^u), where is the gauge field strength 

F% = d,A\ - d v A% - ge\ c A\A\ . (II. 1) 

(For a generic gauge group, e ahc should be replaced with the structure constant of that group.) 
Under the action of gauge transformation U = exp(—X a T a ), where T a are generators of the 
su(2) algebra, 

[T a , T b ] = ie ab c T c , (II.2) 
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A a „ transforms as 

A M -> C/A^t/t _ -[/^^ . (II.3) 

Therefore, out of 12 components of A a ^, nine are physical and three are gauge freedoms, 
which may be removed by a suitable choice of gauge parameter A . Since we are interested 
in isotropic and homogeneous FRW cosmology, the temporal gauge 

A a = 0, (II.4) 

appears to be a suitable gauge fixing. This fixes the gauge symmetry f ill . 3 [) . up to the global, 
time independent SU(2) gauge transformations. This global SU(2) is the key to restoring 
the rotation symmetry in the presence of the background gauge fields. We identify this 
SU (2) with the three-dimensional rotations of the FRW background and, since the physical 
observables of the gauge fields are defined up to gauge transformations (or in other words, 
only gauge-invariant combinations are physical observables) the rotation symmetry may be 
preserved. This latter is done by turning on a specific gauge field configuration in which 
this identification can be made. 

In order to see the above in a more technical language, consider the background FRW 
metric 

ds 2 = -dt 2 + a 2 (t)5 ij dx i dx j , (II.5) 

where i,j = 1,2,3 denote the indices along the spacelike three-dimensional hypersurface 
S, whose metric is chosen to be a 2 Sij. By choosing the (comoving cosmic) time direction, 
metric on £ is then defined up to 3D foliation preserving diffeomorphisms. If we denote the 
metric on constant time hypersurfaces £ by g^, we can introduce a set of three vector fields, 
{e a i (p)}, the triads, spanning the local Euclidian tangent space TS P to the £ at the given 
point p. The triads satisfy the following orthonormality relations 

q ij = e a i e b j 6 ab , 5 ab = e\e h ^ . (II.6) 

The triads are then defined up to local 3D translations and rotations, which act on the 
"local indices" a, b. In particular the triads which are related to each other by local rotations 
A a b G SO (3) 

e\ e\ = A a b (p)e\, (11.7) 

at each point p, lead to the same metric g^. The elements of this local £0(3) may be 
expressed in terms of su(2) generators T a as A = exp(— X a T a ). There are (infinitely) many 
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possibilities for e a i for the FRW metric (1II.5j) . and one obvious choice is 

e\ = a(t)5?, (11.8) 

which identifies the space coordinate indices i, j, k, ■ ■ ■ , with the local frame indices a, b, c, ■ ■ ■ . 

We may now readily identify the remaining global SU(2) gauge symmetry with the global 
part of the 3D rotation symmetry (JTT7TJ) . This can be done through the following ansatz 

a\ = = amm, (n.9) 

where under both of 3D diffeomorphisms and gauge transformations, ip(t) acts as a genuine 
scaler field. Technically, the ansatz (1II.9|) identifies the combination of the gauge fields for 
which the rotation symmetry violation caused by turning on vector (gauge) fields in the 
background is compensated for (or undone by) the gauge transformations, leaving us with 
rotationally invariant background. 

As a result of this identification the energy-momentum tensor produced by the gauge 
field configuration ( 1II.9j) takes the form of a standard homogeneous, isotropic perfect fluid 

T» = diag{-p,P,P,P). (11.10) 

To see this, consider a general gauge and Lorentz-invariant gauge field Lagrangian density 

L = L(F® U ; g^v). The corresponding energy- momentum tensor is given by 

-2 6{y/=gL) _ 5L 
8g^ 5F°, 

To compute T^, we need to first calculate the field strength F^ u for A a ^ in the temporal 
gauge A a = 0, and for the field configuration of flll.9j) : 

(11.12) 

F% = -g^ 2 e% , 

where dot denotes derivative with respect to the comoving time t and for the ease of notation 
we have introduced 

(p = a(t)ip(t) . (11.13) 

(Note that 0, unlike ip, is not a scalar.) It is now straightforward to calculate energy density 
p and pressure P, in terms of <f> and its time derivatives. Plugging All. 12[) into fill. lip yields 

" = &r* _ • £rai - (IU4) 

P = *£4 (11,5) 
7 



T, v = — lV r „,I ' = 27=^^ + g^L ■ (11.11) 



where L re d. is the reduced Lagrangian density, which is obtained from calculating L(F^ V ; g^ v ) 
for field strengths F® u given in (III. 121) and FRW metric (III. 51) . 

One can check that £ rea . is the true reduced Lagrangian for the reduced phase space of 
the field configurations in the ansatz (III.9p (and in the temporal gauge). In order to do this, 
one can show that the gauge field equations of motion 

D *wir = > ( IL16 ) 

where is the gauge covariant derivative, (%) allow for a solution of the form (III.9j) and, (it) 
once evaluated on the ansatz (1II.9P become equivalent to the equation of motion obtained 
from the reduced Lagrangian L red (cf), 0; a(t)) 

, a3 dLr ?L) _dL^ = o ^ (ni7) 



a 3 dt 



In technical terms, there exists a consistent truncation of the gauge field theory to the sector 
specified by the scalar field ip (or 0). In the next section we will study the cosmology of this 
reduced Lagrangian, with a specific choice for the gauge field theory action. 



III. A SPECIFIC GAUGE-FLATION MODEL, ANALYTIC TREATMENT 

In the previous section we showed how homogeneity and isotropy can be preserved in a 
specific sector of any non-Abelian gauge field theory. In this section we couple the gauge 
theory to gravity and search for gauge field theories which can lead to a successful inflationary 
background. The first obvious choice is Yang-Mills action minimally coupled to Einstein 
gravity. This will not lead to an inflating system because, as a result of scaling invariance 
of Yang-Mills action, one immediately obtains P = p/3 and p > 0, and in order to have 
inflation we should have p + 3P < 0. So, we need to consider modifications to Yang-Mills. 

As will become clear momentarily, one such appropriate choice involving F A terms is 



S = d 4 Xy/—g 



R 1 K 

fa p iiv j / ixi/ Act pa pa \2 

2 ^ \xu a ~ gg^ \ ^ u Actv 



(III.l) 



where we have set SnG = M^ 2 = 1 and ^ vXa is the totally antisymmetric tensor. We 
stress that this specific F A term is chosen only for inflationary model building purposes and, 
since the contribution of this term to the energy-momentum tensor has the equation of state 
P = —p, it is perfect for driving inflationary dynamics. The justification of this term within 
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a more rigorous quantum gauge field theory setting will be briefly discussed in section VII. 
(To respect the weak energy condition for the F A term, we choose k to be positive.) 

The reduced (effective) Lagrangian is obtained from evaluating ( IHI.ip for the ansatz (III.9p 



L - d= 2^- — + K ^^ ) - (IIL2) 



The energy density p and pressure P are 



P= 2 { ^ + — + K ^^ ) > (IIL3) 



P=-/- + ^-3^). (III.4) 



As we see p and P have Yang-Mills parts and the F A parts, the k terms. If we denote the 
Yang-Mills contribution to p by p YM and the F 4 contribution by p K , i.e. , 



then 

1 



P = Pym + Pk , P = T^Pym ~ P* ■ ( IIL6 ) 



Field equations, the Friedmann equations and equation of motion, are then obtained as 

H = n ( + — T" + ~r ) > ( IIL7 ) 



* = + (HIS) 

(1 + /t— — ^— ) — + (1 + «5L ^f - + (1 - 3^)^ = . III.9 
ar a a z a 6 ar a 

We start our analysis by exploring the possibility of slow-roll dynamics. To this end it is 
useful to introduce slow-roll parameters 1 

EH , 

e = ~lp- <^-WJi' (IIU0) 

where e is the standard slow-roll parameter and r) is related to the time derivative of e as 

" = f -drT (mil) 



1 We note that our definition of slow-roll parameters e, 77 for the standard single scalar inflationary theory 
L = y 2 - V{ V ) reduces to fl e = (v) > »? = M £tt- 



Therefore, to have a sensible slow-roll dynamics one should demand e, rj <C 1. Using the 
Friedmann equations (IIII.7j) and (IIII.8I) and definitions flIII.5j) we have 

e = 2pYM . (111.12) 

Pym + P« 

That is, to have slow-roll the K-term contribution p K should dominate over the Yang-Mills 
contributions p YM , or p K ^> p YM . As we will see, the time evolution will then increase 
p YM with respect to p K , and when p YM ~ p K , the slow-roll inflation ends. Noting that 
p + 3P = 2(p YM — p K ), inflation (accelerated expansion phase) will end when p YM > p K . 

For having slow- roll inflation, however, it is not enough to make sure e < 1. For the 
latter, time-variations of e and all the other physical dynamical variables of the problem, 
like rj and the if) field, must also remain small over a reasonably large period in time (to 
result in enough number of e- folds). To measure this latter we define 



5 = 



Htp 



in terms of which the equations ( 1IH.7l) -( lHI.9l) take the form 

e = 2 

77 = e - (2 - e 



5 



HI 



35 

5)~e + T 



(111.13) 

(III. 14) 
(111.15) 



Comparing (lIH.llj) and (1IH.15|) we learn that to have a successful slow-roll, e ~ He 2 and 
rj ~ e, we should demand that 5 ~ e 2 . Explicitly, the equations of motion (IIII.7j) . flIII.8j) and 
f ]III.9j) admit the solution 2 

e ~ V ;2 (7 + 1)> 



V 



^ (2-6)( 7 + l) 3 

2 S ' 

9 £ 6 

where ~ means equality to first order in slow-roll parameter e and 3 



7 



i/ 2 



or equivalently H 



9 2 e 



e-ip 2 7(7 + 1) 



(111.16) 
(111.17) 

(111.18) 



(111.19) 



2 Our numerical analysis reveals that even if we start with S/(H5) ~ 0(1), while ipf ~ e <C 1, after a 
short time it becomes very small and hence for almost all the inflationary period we may confidently use 
e 2 . See section HV] for a more detailed discussion. 



u - 6( 7 +l) C 

3 Note that all the dimensionful parameters, i.e. H,ip and k, are measured in units of M p \; if, tp have 
dimension of energy while k has dimension of one- over- energy density. 
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In the above 7 is a positive parameter which is slowly varying during slow-roll inflation. 
Recalling 13j) and that 5 ~ e 2 , flT 



I.19j) implies that 7 if 2 remains almost a constant 



during the slow- roll inflation and hence 1| 

e _ 7 + 1 7 _ #? 



(111.20) 



where e^, 7$ and ifj are the values of these parameters at the beginning of inflation. As 
discussed the (slow-roll) inflation ends when e = 1, where 

7< + 1 ^/ T< /ttt 01 \ 

7/ - —— , 77f - — T T e i ■ ( IIL21 ) 
Q ti i 7i + 1 

Using the above and (IIII.10I) one can compute the number of e-folds N e 



IV. NUMERICAL ANALYSIS 

As pointed out, our gauge-flation model has two parameters, the gauge coupling g and 
the coefficient of the F A term k. The degrees of freedom in the scalar sector of the model 
consists of the scalar field ip and the scale factor a(t) and hence our solutions are specified by 
four initial values for these parameters and their time derivatives. These were parameterized 
by Hi, €i and ipi and 5i (or ipi). The Friedmann equations, however, provide some relations 
between these parameters; assuming slow-roll dynamics these relations are (1III.16I) - (1III.18I) . 
As a result each inflationary trajectory may be specified by the values of four parameters, 
(g, K] ipi, ipi). In what follows we present the results of the numerical analysis of the equations 
of motion ( 1III.7I) . ( 1III.8I) and ( 1III.9I) . for three sets of values for (^,^5 g, k). 



► Discussion on diagrams in Fig[li 

The top left figure shows evolution of the effective inflaton field ip as a function of Hit. 
As we see, there is a period of slow-roll, where ip remains almost constant and e is almost 
constant and very small. Toward the end of the slow-roll e grows and becomes one (the 
top right figure), (slow- roll) inflation ends and ip suddenly falls off and starts oscillating. 
As we see from the top right figure, the slow-roll parameter e has an upper limit which is 
equal to 2. This is understandable recalling (1111.121) and that p K is positive definite. At 
the end of slow- roll inflation p K , is negligible and the system is essentially governed by the 
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Yang-Mills part p YM - In addition, the top left figure shows that amplitude of the tp field in 
the oscillatory part is dropping like t~ 1//2 . (The minima of e is also fit by a t l l 2 curve.) This 




FIG. 1: The classical trajectory for ipi = 0.035,^ = -1(T 10 ; g = 2.5 x 1(T 3 ,k = 1.733 x 10 14 . 
These values correspond to a slow-roll trajectory with Hi = 3.4 x 10 -5 , = 6.62, ej = 
9.3 x 10~ 3 , 5i = 8.4 x 10~ 5 . These are the values very close to the range for which the gauge- 
flation is compatible with the current cosmological and CMB data (c/. discussions of section |VT|) . 
Note that k, Hi and ipi are given in the units of M p \. 

behavior is of course expected, noting (IIII.5I) and that in the oscillatory regime the dominant 
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term is the p YM , that is, the system effectively behaves as a g 2 ip 4 chaotic inflation theory. 
And it is well-known that after the slow-roll phase ip(t) in the g 2 ip A theory oscillates as a 
Jacobi-cosine function whose amplitude drops like t~ x l 2 13|. In other words, the averaged 
value of e and a(t) behave like a radiation dominated Universe (recall that for a radiation 
dominated cosmology e = —H/H 2 = 2). 

The bottom left figure shows the phase diagram of the effective inflaton trajectory Note 
that this diagram depicts 4>/a(t) vs <f)/a(t) (rather than ip vs if}). The rightmost vertical line 
is where we have slow-roll, because = a(t)ip and during slow-roll if) is almost a constant. 
The curled up part is when inflation has ended, and when the system oscillates around a 
radiation dominated phase. This latter may be seen in the figure noting that the amplitudes 
of oscillations of both <j>/a(t) and (f>/a(t) drop by t~ x l 2 . The bottom right figure shows number 
of e-folds as a function of comoving time. As expected, the number of e-folds reaches its 
asymptotic value when e ~ 1. 

One can readily check that the behavior of if), e, and the number of e-folds during slow-roll 
inflationary period has a perfect matching with our analytic results of previous section. We 
note that, as will be discussed in section fVTl the set of parameters Hi, ji, if>i, g corresponds to 
an inflationary model close to the range of values compatible with the current cosmological 
and CMB data. 

► Discussions on diagrams in Figs. \M and0 

Fig. [2] corresponds to a slow-roll trajectory which starts with a lower value of e, but 
almost the same value for H, compared to the case of Fig. [TJ For this case we hence get 
a larger number of e-folds. The qualitative shape of all four figures is essentially the same 
as those of Fig. HJ and both are compatible with our analytic slow-roll results of previous 
section. Our numeric analysis indicates that the behavior of the phase diagram for if> field 
and e do not change dramatically when the orders of magnitude of the initial parameters 
are within the range given in Figs. [TJ or [2] 

Fig. [3] shows a trajectory with a relatively large if). As we see after a single fast falloff the 
field falls into usual slow-roll tracks, similar to what we see in Figs. [Hand [2J The oscillatory 
behavior after the inflationary phase, too, is the same as those of slow-roll inflation. The 
graph in the square in the bottom left figure shows, with a higher resolution, the upper part 
of the phase diagram which comes with under brace. This part corresponds to the dynamics 
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FIG. 2: The classical trajectory for tpi = 0.025,^ = -HT 10 ; g = 2.507 x 1CT 3 ,k = 1.3 x 10 15 . 
These values correspond to a slow-roll trajectory with Hi = 3.63 x 10~ 5 , = 2.98, e, t = 
2.5 x 10~ 3 , 5i = 1.1 x 10~ 4 . These figures show that it is possible to get arbitrarily large numbers 
of e-folds within the slow-roll phase of our gauge-nation model. 

of i/j field after inflation ends, and as expected has the same qualitative form as the phase 
diagrams in the slow-roll trajectories of Figs. [1] and |2j Our numeric analysis shows that 
getting a large enough number of e-folds and the generic after-inflation behavior of the fields 
is robust and does not crucially depend on the initial value of the fields, ^ or but it is 
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FIG. 3: The classical trajectory for tpi = 8.0 x KT 2 ,-^ = -1CT 4 ; g = 4.004 x 10~ 4 , n = 4.73 x 10 12 
These values correspond to a non- slow-roll trajectory with 5 ~ 2, Hi = 6.25 x 10 -4 , £j = 6.4 x 10 -3 . 
We start far from the slow-roll regime for which 5 ~ e 2 1. This latter is also seen from the 
phase diagram (bottom left figure). Despite starting far from slow-roll regime, as we see from 
the top left figure, after an abrupt oscillation the field ip loses its momentum and falls into the 
standard slow-roll trajectory. As shown in the bottom right figure, for this case we get a large 
number of e-folds. Getting a large enough number of e-folds seems to be a fairly robust result not 
depending much on the initial value of 5. 



15 



sensitive to the initial value of e. More precisely, as long as e remains small of order 0.01, 
regardless of the value of 5 we can get arbitrarily large number of e-foldings. The examples 
of small and large S values have been, respectively, given in Figs. |2] and El 

It is also useful to work out the displacement of the ip field during inflation. To this end, 
let us start from ip = —5Hip and use the value for 5 given in (JTTT7T7J) . Our numerical analysis 
reveals that the dynamics of the system is such that even if we start with of order 1 — 10, 
but ipf ~ €i, after a short time Hit ~ 1, 5/(H5) becomes very small and hence in almost all 
the inflationary period, except for the first one or two e- folds, 5 ~ 6 ^ +1 ^ 2 - Since variations 
of 5 in this period happens very fast, and after that it remains almost zero, during this 
period ip does not change much. Numerical analysis also shows that these arguments are 
generically true even if we start with a large value of 5, as in Fig. [3J provided that the other 
parameters are such that we get large number of e-folds (about N e ~ 60 or larger). This 
latter can be easily arranged for. This is again confirming the robustness of the classical 
inflationary trajectories with respect to the choice of the initial conditions. Therefore, we 
may confidently compute roaming of ip field using the equation 

ip ^ g 2 H 
~ ~6~FP ' 

Integrating the above equation one can compute the displacement of the scalar field ip during 
inflation. If we denote the value of ip in the beginning and end of inflation, respectively, by 
ipi and ipf we obtain 

^?-^ 2 . (iv.i) 

Alternatively one could have used ( 1111.181) ng 2 ip 6 (l — 5) 2 = 2 — e to compute the change in 
ip. If we are in slow- roll regime where we can drop 5-term, then the ratio of ijjj (which is 
computed for ej = 1) to ipi is obtained as ^ = And up to percent level, y/3/2 and 

are equal, a confirmation of the validity of the slow-roll approximations we have used. 
Interestingly, at this level of approximation, the roaming of the ip field is independent of the 
initial value of the 7 (or the initial value of Hubble) parameter and ipi and ipf have the same 
orders of magnitude. 
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V. COSMIC PERTURBATION THEORY IN GAUGE-FLATION 



So far we have shown that our gauge-nation model can produce a fairly standard slow- 
roll inflating Universe with enough number of e-folds. The main test of any inflationary 
model, however, appears in the imprints inflation has left on the CMB data, i.e. , the 
power spectrum of curvature perturbations and primordial gravity waves, and the spectral 
tilt of these spectra. To this end, we should go beyond the homogeneous (^-independent) 
background fields and consider fluctuations around the background. This is what we will 
carry out in this section. 



A. Gauge- invariant metric and gauge field perturbations 

Although not turned on in the background, all of the components of metric and the 
gauge field in all gauge and spacetime directions will have quantum fluctuations and should 

2 



,13] 
(V.l) 



be considered. The metric perturbations may be parameterized in the standard form 

ds 2 = -(1 + 2A)dt 2 + 2a(diB - SJdx'dt 

+ a 2 ((1 - 2C)5 ij + 2d tj E + 2d (i W j ) + dx l dx j , 

where di denotes partial derivative respect to x l and A, B, C and E are scalar perturbations, 
Si, Wi parameterize vector perturbations (these are divergence-free three- vectors) and hy, 
which is symmetric, traceless and divergence-free, is the tensor mode. The 12 components 
of the gauge field fluctuations may be parameterized as 

5A\ = 5%Y + ^ , (V.2) 
5A\ = 8?Q + 6 ak d tk M + g<\>e\ k d k P + 8id iVj + e°>,- + 5 a % , (V.3) 

where, as discussed in section [Til we have identified the gauge indices with the local Lorentz 
indices and the expansion is done around the background in Aq = temporal gauge. In 
the above, as has been made explicit, there are four scalar perturbations, Q, Y, M and P, 
three divergence-free three- vectors Ui, i>, and Wi, and a symmetric traceless divergence-free 
tensor Uj, adding up to 4 + 3 x 2 + 2 = 12. Q is the perturbation of the background field 
<p, which is the only scalar in the perturbed gauge field without spacial derivative. We are 
hence dealing with a situation similar to the multifield inflationary theories where we have 
adiabatic and isocurvature perturbations. If the analogy held, Q would have then be like 

17 



the adiabatic mode. However, as we will see this is not true and the curvature perturbations 
are dominated by other scalars and not Q. As another peculiar and specific feature of the 
gauge-flation cosmic perturbation theory, not shared by any other scalar-driven inflationary 
model, we note that the gauge field fluctuations contain a tensor mode £y. As we will show, 
the power spectrum of this mode is nonzero, but its effect is such that the power spectrum 
of the tensor modes ends up to be exactly those of usual scalar-driven inflationary models, 
with only the hij of the metric contributing to the tensor power spectrum. 

Because of the presence of gauge symmetries not all 10 + 12 metric plus gauge field per- 
turbations are physical. Altogether there are four diffeomorphisms and three local gauge 
symmetries, hence we have 15 physical degrees of freedom. The four diffeomorphisms re- 
move two scalars and a divergence-free vector [2J, and the three gauge transformations one 
scalar and one divergence-free vector. Therefore, we have five physical scalar perturba- 
tions, three physical divergence-free vector, and two physical tensor perturbations. (These 
amount to 5x1 + 3x2 + 2x2 = 15 physical degrees of freedom.) The gauge degrees 
of freedom may be removed by gauge-fixing (working in a specific gauge) or working with 
gauge-invariant combinations of the perturbations. In what follows we work out the gauge- 
and diffeomorphism-invariant combinations of these modes. 

1. Scalar modes 

Let us first focus on the scalar perturbations A, B, C, E, Q, Y, M and P. Under 
infinitesimal scalar coordinate transformations 



where St determines the time slicing and Sx the spatial threading, the scalar fluctuations of 
the gauge field and metric transform as 



t-tt = t + 5t, 



(V.4) 



x l ->■ x l = x l + S l3 djSx , 



Q->Q-<f>5t, 



Y->Y-<f>5x, 



M -> M - 4>Sx , 



P -> P, 



A->A-6t, 
B ->■ B + — - aSx , 



C ^C + HSt, 



(V.5) 



E ->■ E-Sx. 



a 
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On the other hand under an infinitesimal gauge transformation A a , fluctuations of the 
gauge field transform as 

5A\ 5A\ - U,\ a - e\ c \ b A c ^ . (V.6) 

The gauge parameters A a can be decomposed into a scalar and a divergence-free vector: 

A a = 5 ai di\ + 6?\j . (V.7) 

The scalar part of the gauge field perturbations under the action of the scalar gauge trans- 
formation A transform as 

Q^Q, Y^Y--\, 

1 i (V.8) 

M > M • -A. P^P+-\. 

9 9 

We note that Q is gauge-invariant and this is a result of identifying the gauge indices with 
the local Lorentz indices and that Q is a scalar. 

Equipped with the above, one may construct five independent gauge-invariant combina- 
tions. One such choice is 4 

= C + a 2 H(E--), (V.9) 

CI 

• = *-s ("*<*-!>)• (V10) 

Q = Q - a 2 <j)(E - — ) , (V.ll) 

CI 

M = M + P-(f)E, (V.12) 
y = Y + P-(f)E. (V.13) 

The first two, $ and \P are the standard Bardeen potentials, while Q, M and V are the three 
gauge and diffeomorphism-invariant combinations coming from the gauge field fluctuations. 

Finally, for the later use we also present the first-order perturbations of the gauge field 
strength sourced by the scalar perturbations 

5E a 0l = 5?Q + 5 aj dij(M -Y)- ge^d^P) + <f>Y) , (V.14) 
6F% = 25 a [] d l] {Q + g 2 < p 2 P) + 2g < Pe%d, ]k {P + M) -2<?e%.Q0. (V.15) 

Note that are not gauge- invariant, as under gauge transformations F° — > F^ v — 



4 These choices are not unique and one can construct other gauge invariant combinations too. 
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2. Vector modes 



Next, we consider the vector modes Si, Wi, Ui, v% and Wi. Under infinitesimal "vector" 
coordinate transformations 

x* -> x* = x* + 6x' v , (V.16) 

where diSxf = 0, 

Si -»■ Si + adxl , Wi^rWi- 8x1 , 

( V - 17 ) 

Ui —7- Ui — (j)0X v , Vi — > Vi — (j)SXy , Wi — > Wi . 

On the other hand under the vector part of infinitesimal gauge transformation (IV. 61) . 

Ui ^ Ui - -\y , Vi -> Vi --Ay , Wi-^Wj + ^Ay, (V.18) 

and obviously Wi remain invariant. 

The three gauge- and diffeomorphism-invariant divergence-free vector perturbations may 
be identified as 

Z t = aWi + Si, (V.19) 
U, = Ul - Vl + 4>Wi , (V.20) 
V, = Wl + - 0Wi) . (V.21) 

The contribution of vector perturbations to the first-order gauge field strength perturba- 
tions are 

=Kdi(vj - Uj) + eV (wj + 94>Uj) , 
<fF%. =2e a y + «^ fc ) + 2#5 a [iWj] . 

3. Tensor modes 

One can show that the tensor perturbations hij and ty, being symmetric, traceless and 
divergence- free, are both gauge- and diffeomorphism-invariant. The contribution of to 
the first order perturbed corresponding to is 

^ = ' % ' (V.23) 
5F£=<5 ofe %t i]fc -#e^ ]fc . 
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B. Field equations 



Having worked out the gauge invariant combinations of the field perturbations, we are now 
ready to study their dynamics. These first order perturbations are governed by perturbed 
Einstein and gauge field equations 

5G, U = 5T, U1 5 (d^-^=0, (V.24) 

where by 5 in the above we mean first order in field perturbations. Since we are dealing 
with an isotropic perfect fluid in the background, as it is customary in standard cosmology 
text books j^, it is useful to decompose energy- momentum perturbations as 

ST i:j =P 6 gij + a 2 (<%£P + dijn 3 + durj + djirf + vrg) , (V.25a) 
5T l0 =P 5g l0 - (P + Po )(di5u + 8uJ) , (V.25b) 
ST 00 = - p 5g 00 + 5p , (V.25c) 

where subscript, "0" denotes a background quantity and tt s , ttY , vr^ represent the 
anisotropic inertia and characterize departures from the perfect fluid form of the energy- 
momentum tensor; tt^ and nfj and the vorticity 5uf satisfy 

di5uY = d&Y = 0, vrj = 0, 9^5 = 0. (V.26) 

Since being a perfect fluid or having irrotational flows are physical properties, their corre- 
sponding conditions are gauge-invariant. 



1. Scalar modes 

As is usually done in cosmic perturbation theory, it is useful to write down the equations 
of motion in a gauge-invariant form. In order this we note that ST^ has four gauge-invariant 
scalar parts 5p g , 5P g , 5q g , 

6p g = 6p-p a 2 (E-^), (V.27) 
5P g = 5P- P a 2 {E - -) , (V.28) 
5q g = 5q+ (p + P )a 2 (E - -) , (V.29) 
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and Ti s [2J, where Sq = (po + Po)Su. Out of the 10 perturbed Einstein equations, there are 
four scalars, two (divergence-free) vectors, and one massless tensor mode (gravitons). Among 
the four scalar perturbed Einstein equations, one is dynamical and three are constraints. 
Therefore, they do not suffice to deal with five gauge- invariant scalar degrees of freedom and 
one equation is missing. This last equation is, of course, provided by the perturbed gauge 
field equations of motion. To this end, we use two observations: 

• Using ( 1V.14P and after a lengthy algebra (which was performed by Maple codes) one 
can show that the anisotropic inertia n s is given by the following linear combination 
of the gauge-invariant quantities 5 

a ' a 3 
We can write ^ in terms of a 2 ir s and the rest of variables. 



a V = 2^(aV - aM) + 2^aM . (V.30) 



Computing the second-order action for the gauge field perturbations, we observe that 
the momentum conjugate to the new variable a 2 ir s , is identically zero and hence a 2 ir s 
is a nondynamical variable; it is a constant of motion. To see the latter, note that 
there is no V in the SF a (cf. (jV.14p ). In our analysis we can hence consistently 



choose the initial conditions so that 7r s vanishes. 

The above will then provide the last missing equation, and the perturbed Einstein equations 
are enough to solve for the scalar perturbations. The other perturbed gauge field equations 
of motion do not lead to independent equations. Studying a generic constant a 2 ir s is an 
interesting question we postpone to future works. 

nn 

With the above argument we have [2|, |3J 

aV = $ = (V.31) 



5 Note that in the first order perturbation theory the tensor and vector modes do not contribute to the 
scalar parts of the energy-momentum perturbations. 
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We hence remain with three equations for the three variables. Using (1V.14[) . and after 
lengthy calculations (confirmed by Maple codes too), we obtain 



5 Ps 



3(1 



Kg 



2^4 



(j)Q k 2 



a 4 ' a 2 



Wl K<P\g 2 (t)\ T Q v 



5P n 



:i-3 



Kg 



2^4 



a 1 



a 2 3a 2 



(1 -3^ ) + 2(1 -3^ 



,2^,3 



■oM 



,2 A3 



■oM 



2(l-3^)^(* + ?) 



a 2 ' a 4 



0Q ^ 2 3 
a 2 
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4^6 



a 3 a 5 



(V.32) 



(V.33) 
(V.34) 



Note that in the above expressions tt s = condition and f ]V.31j) have been employed. 

We are now ready to write down the three perturbed Einstein equations, two of which 
are constraints and one is dynamical: 



6q g + 2(tf + HV) = , 

k 2 

5p g - 3H5q g + 2—^ = , 

5P g + 5q g + 3H5q g + (p + P )* = • 



(V.35) 
(V.36) 
(V.37) 



In the above we have already used $ = ^ relation. These relations provide enough number 
of equations for the gauge-invariant scalar perturbations to which we return in the next 
subsection. 
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2. Vector modes 



To study the vector perturbations, we first work out vector parts of the perturbed energy- 
momentum tensor, bq( and , using f)V.22[) : 6 

^2 



8ql = u>i H (aui + -aS t ) + 

a z \ a a I a z 



i;2 



+ ^(Vx (j-u)) .-4(Vx {vS + g^t , (V.38) 

Ct \ /id 

7r/ = — 5~( 1 Vi -Wi) + -{u i -v i + -W i ), (V.39) 

a a a a a a 

which as expected, are gauge-invariant and hence can be written in terms of physical gauge 
and diffeomorphism-invariant variables as 



6q? = (-Vi Vi H -allj H T aZi) - V x ( 1 -U) , (V.40) 

a a a a z a A \ a a a 1 



TTj = —5 h -Ui . (V41) 

a z a a 

The perturbed Einstein equations have two vector equations, one constraint and one dy- 
namical equaiton. These equations are 

di f 2a 2 7r/ - = , (V.42) 

k 2 

2SqY + —aZi = 0. (V.43) 
a z 

To solve for the three gauge-invariant vector perturbations the two Einstein equations 
are not enough and we should use information from the perturbed gauge field equations. 
As for the latter, we note that one can use (IV. 391) to write Ui in terms of a 2 nY and the 
other variables. On the other hand, since there is no iii in (IV.22[) . once we write down 



the second-order action for the gauge field perturbations the momentum conjugate to tiJ is 
vanishing and hence this variable is nondynamical; it is a constant of motion. We may then 
choose the initial conditions such that 

a 2 vrf = . (V.44) 



6 Note that at first order in perturbation theory only the vector perturbations contribute to the vector part 
of energy- momentum tensor perturbations. 
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This completes the set of equations we need for solving vector perturbations. Eq. (1V.42j) 
then implies that 

Zi = 4 ■ (V.45) 

The above is the usual result of the scalar-driven inflationary models that the vector modes 
are diluted away by the (exponential) accelerated expansion of the Universe during inflation. 
In our model, despite having vector gauge fields as inflaton, we confirm the same result. 7 



3. Tensor modes 



As discussed, there are two gauge- and diffeomorphism-invariant tensor modes and 
tij, while perturbed Einstein equations only lead to one equation for h^. This equation, 
which is sourced by the contribution of to the energy-momentum tensor, reads as 



+ 3Hhij + —hij = 2%^ 



(V.46) 



The other equation of motion is provided with the perturbed gauge field equations of motion. 
After a tedious but straightforward calculation, which is also confirmed by the Maple codes, 
we obtain the following second-order action for the tensor modes 



5S^ ] ~ - / d 6 xdt a 



2 { a 2 



ar J a a a 6 a 



1 k 2 

a? K v a 2 tJ 



K,g 2 (f) 2 4> 2 2 l- 2 k 2 2 

; 2 2 ij' ~A\ ij 2 ijt 

a z a z J 4 J cr J 



Note that in the above we have already used the slow-roll approximation 



(V.47) 



H(f). From 



the above second order action one can readily compute ir. 



■ij 



a 4 a a a A a 



(V.48) 



7 In contrast to standard scalar driven inflationary models, in gauge-flation there is the possibility of a 
non-zero but constant a 2 irY which leads to 



2aX 
aH 



c 



The above, still shows a suppression by a 1/a factor (to be constrasted with standard 1 /a 2 suppression) 
and hence again the power spectrum of the vector modes are unimportant in inflationary cosmology. 
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Being traceless and divergence-free Uj and each has 2 degrees of freedom which are 
usually decomposed into plus and cross (+ and x) polarization states with the polarization 
tensors e^' x . Since we have no parity- violating interaction terms in the action the equations 
for both of these polarization have the same time evolution and one may then introduce h 
and T variables instead 



and in the computation of the power spectrum consider these variables, treating them as 
scalars, but multiplying by a factor of 2 to account for the two polarizations. The second- 
order action for h and T reads as 



SSP = l - I <f\nlr 



\^ 2l K 2 {l - x)h 2 + 2^ 2t K 2 {~\- + xT)h 
2 Jx 



+ ^ 2 {T' 2 - k 2 T 2 - X 2 e^J-T 2 ) + - (h' 2 + (2 - e)jK 2 h 2 - k 2 h 2 ) 



(V.50) 



where r is the conformal time dt = adr, < K = a and prime denotes a derivative with respect 
to the conformal time. 

C. Primordial power spectra and the spectral indices 

In the previous part, we provided the complete set of equations which govern the dynamics 
of scalar, vector and tensor modes. In this subsection we set about solving these equations, 
quantize their solutions, and compute the power spectra. 

1. Scalar modes 

In order to determine the power spectrum of the scalar perturbations we have to deal 
with two constraint (IV.35j) . flV.36j) and one dynamical equation flV.37j) . In contrast to the 
case of single scalar field inflationary models which one of the constraints (a combination 
of 5P and Sq equations) reduces to the equation of motion of the background field, in our 
case both of them remain independent and should be considered. In Appendix C we provide 
more details about this issue. The constraint equations (lV.35j) and (lV.36j) in the slow-roll 
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regime take the form 

IK U> - il>~ + ^V 3 (7 + l)aMJ + (* - ^ 3 (aM))' = , (V.51) 
IK I — + (1 + — ^-)a* + — (1 + — ;r-)a* -(! + - + ^ 2 (1 - -))<?V 



V> 2 J 3 V 2 7 7 ^ 

+IK 2 (2 + ^ 2 ( 7 - l))(a* + -) = . (V.52) 

where ' denotes derivative with respect to the conformal time r. 

To analyze the dynamical equation (TV.37j) . we note that for a general hydrodynamical 
fluid with pressure P(p, S), the pressure perturbation 8P can be decomposed as 

5P = c 2 s 8p + 78S , (V.53) 

where 5S is the entropy perturbations, c 2 = (|p)s is the speed of sound, and 7 = (§f V 
With this decomposition, combining ( IV. 36|) and ( IV. 37|) yields 



+ 3IK(1 + c 2 )^' + c 2 s k 2 m + (2IK' + (1 + 3c 2 )IK 2 ) = ^a 2 T55 , (V.54) 



which upon introducing 



9 l/T f« ■■ * (V.55) 
where po and Po are background energy density and pressure, simplifies to 

nil 

u" + c 2 A; 2 m - — u = 8s , (V.56) 

where c 2 = Next, we need to compute 8s explicitly. We will do this in two regions, the 
asymptotic past kr — > — oo and the supperhorizon region kr — > 0. 

► The asymptotic past behavior of 8s. 

To study the effects of 8s term in the asymptotic past kr — > — oo (Minkowski) limit, we 
note that in this limit constraints (IV. 511) and (IV. 521) take the following forms 

Q-^(|-^ 2 aM)' = 0, (V.57) 

Q' - ^-(1 + -)g 2 i, 2 aM + ^ ( ^ - (1 - -)g^ 2 aM) = . (V.58) 

We solve (IV. 58|) by considering the following ansatz 8 

V - ^(1 + -)9 2 ^ 2 aM = ^ 2 Q', (V.59) 
oil 7 Z 



The validity of the ansatz is proved in Appendix B. 
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with W being a constant number of order one, to be determined momentarily. In this limit 
from (IV.32j) . (IV.33j) and flV.53j) one can read the value of 5 s 



5s 



4(- )Hip- 



Q' 



From flV.591) . the constraint equation (IV.58j) takes the form 

ry 1 9k 2 

(W - 2- )HibQ' + ^ ~ , 

v 7 + 2 ; v 3 



(V.60) 



(V.61) 



and hence 



5s 



V 7-1 

3 ( 7 + 2)W - 2(7 - 1) 

( 7 + 2)W+ |(7-1) 



u 



3 " 2( 7 - 1) - ( 7 + 2)W 
On the other hand, from (IV. 571) . (1V.59I) and (1V.61j) we learn that 

2 2 7 + 2 



(V.62) 
(V.63) 



3 3W( 7 + 2) + 2' 

Consequently from equating the above two values for the speed of sound, we find a quadratic 
equation for W whose solutions denoted by W± are 

2 7-l . 2 



-2 



3 7 + 2 
7 + 1 



,2 



7 + 2 



1 7 + 2 

'3 7 



(V.64) 
(V.65) 



As we see speed of sound squared for one of the solutions is one and for the other one, 
recalling that 7 is positive by definition, is always negative. We also note that all the above 
analysis has been carried out in the slow-roll regime, and in the Minkowski limit kr — > —00. 

► The superhorizon behavior of 5s. 

We now turn to the question of large scale superhorizon behavior of 5s in kr — > 0. In 
this limit, from flV.32j) . (IV.33j) and (IV.53j) one can read the value of 5s 

Q/a , Of /a 



5s = 2^V(27(m + 



+ 



>)■ 



On the other hand, constraint (IV.52I) in the kr — > can be written as below 

(- + a^)' + 2JC(- + a^) =0(e), 



(V.66) 



(V.67) 
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which implies that H + a\l/ consists of a damping term proportional to a 2 and a part of the 



order e. Hence, using ( TV. 6 71) . ignoring the damping terms and up to the first order of e, 5s 
takes the form 



5s = -2K 2 ^ 2 



ft 



it 



(V.68) 



which is a small quantity of the order e. 

Upon using (IV.64I) . (IV.65j) and (1V.68[) . the equation for u (1V.56|) leads to two equations 
for each of the values of speed of sound c 2 ±: 



u'l + (k 2 - CK 2 )u + ~ , 

u »_ _(Il±^ + Cft>_~0, 
3 7 

where, as before, ~ means first order in slow-roll parameters and 

^~ 6 ^ " l 7 + l + He 2 ' 7 + 1 



(V.69) 
(V.70) 

(V.71) 



In the first order in slow-roll parameter e one may ignore the time-dependence of e during 
slow-roll expansion, and 



1 



(1 - e)ft 



and hence 



(ft 2 



K. 



(V.72) 
(V.73) 



To summarize, 5s is of order one in the asymptotic past kr — > — oo limit and is of order 
e in the superhorizon kr — > limit, effecting the spectral tilt (cf. (1V.71j) ). 

The general solutions to ( 1V.69I) and ( IV. 701) for u± can be expressed as a linear combination 
of Hankel and -f/i 2 '', and modified Bessel functions l v and K u . Recalling (1V.55I) . this 
leads to the following solutions for \T/ 



V + (k,T] 



[btHV(k\T\) + b$HV(k\T\)] 



2k 



V^k 



(V.74) 



3 11/' Z *""U \ I / Q 

Before moving onto considering Q equations, some remarks are in order: 
• It can immediately be seen that \T/± given in (IV.74[) . in the asymptotic past limit 



kr — > — oo, indeed reproduce the solutions to 

u" + k 2 c 2 s± u ~ 0. 
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As such, we have found solutions which interpolate between the superhorizon regime 
kr — > and deep subhorizon regime kr — > — oo. 

We stress that as we see 5s in ( IV. 621) and (IV. 681) . in both kr — > — oo and kr — > 
regimes, only depends on u, the variable which plays the role of adiabatic perturba- 
tions. Therefore, unlike the multifield inflationary models, 5s does not represent an 
independent "entropy" mode. In other words, despite having two scalar modes Q and 
M, we do not have entropy perturbations in our system. This is due to the fact that, 
unlike the usual two-field inflationary models, we have two independent constraints 
(IVToTj) and (TV52|) . rather than a single constraint in the usual two-field models jl^j], re- 



lating Q and M. Therefore, we have a single adiabatic perturbation, as in the standard 
single field models. 9 

Eq. ( 1V.74j) is not exact in the sense that the constraint equations, the result of which 
have been plugged into the u equation ( IV. 561) . have only be solved in the kr — >■ and 
kr — > — oo regimes. 

Eq. (IV.69j) implies that u + admits two physical regimes: oscillating regime in the 
asymptotic past kr — > — oo, and the superhorizon limit kr — > where u + freezes out. 
As given by ( IV.70p . m_ has an exponentially damping regime in the asymptotic past 
Minkowski limit kr — > — oo, while it freezes out, just like u + , in the superhorizon 
kr — > regime. The other two scalars, Q and oM, have the same generic behavior. 

Note also that the coefficients bf, i = 1,2 are not completely fixed from the above 
considerations and to determine them we need to know Q, explicitly Q + , to which we 
will turn now. 



Classical solutions for Q. 

To find the equations of motion for Q one may use the solutions ( ]V.74j) and ( 1V.57j) . (1V.58j) . 



Alternatively, one can work out the second-order action for the perturbations Q, \1/ and M 
and insert the constraint equations ( IV. 5 71) and ( 1V.58I) into the action. The second order 



action computation is indeed very tedious, lengthy and cumbersome, but that is necessary for 



9 Recall that the situation here is indeed very much like the single scalar case, where the entropy perturba- 
tions, like ours, is of order slow-roll parameter e. 
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quantization of the perturbations. This is because for performing the canonical quantization 
of the modes, besides the equations of motion we need to have the canonical (conjugate) 
momentum too. In the Appendix A we have presented the explicit form of the second-order 
action, after imposing the gauge-fixing conditions and setting tc s = 0. 

Appearance of negative c 2 modes may cause a concern about a possibility of ghost insta- 
bility in our system. Theoretically we do not expect finding ghosts in our theory because, i) 
we are dealing with a gauge-invariant action and we respect this gauge symmetry. (To be 
more precise, it is spontaneously broken by the choice of classical inflationary background. 
However, as is well-established, spontaneous gauge symmetry breaking does not lead to a 
break-down of Slavnov- Taylor identity which reflects the gauge symmetry and its conse- 
quences about renormalizability and unitarity.) ii) Although we are dealing with a "higher 
derivative" action flHI.lj) . the higher derivative term has a special form: it does not involve 
more than time-derivative squared terms. (This fact is also explicitly seen in (1III.9j) in that 
the equation of motion does not involve more than second time derivative.) As such 
we expect not to see ghosts usually present in the higher derivative theories. Besides the 
above arguments, to make sure about the absence of ghosts, we have explicitly computed 
the second-order action. The expression for the second-order action, after implementing the 
constraints ( 1V.57j) and ( 1V.58|) . explicitly shows that neither Q nor M has negative kinetic 
terms and hence there is no ghost instability in our system. The explicit expression for 
the second-order action is presented in Appendix A and here we only present the simplified 
result for the canonical momenta and for the equations of motion for Q. 

The Q equations of motion in the slow-roll approximation, after using (IV. 591) . (IV. 6411 and 
( 1V.65|) . and some lengthy algebra is obtained to be 

Q'j + (k 2 - (2 + 3C - 4e)J£ 2 ) Q+ ~ , 

fx + 2 \ ( V - 75 ) 

Q" - ( ?^k 2 + (2 + 3(-4e)X 2 ) Q_ ~ , 

where Q± in an obvious notation is related to values of W± and c 2 ± given in ( IV. 64ft and 
(1V.65I) . Defining vq as 

4- 1 = (2 + 3C-4e):KV ~ 2 + 3C, (V.76) 
leading to uq ~ | + (, (1V.75j) takes the form of a standard Bessel equation with the general 
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solutions 



Q + (k,r) ~ V^H e ^)f (*| T |) +dtHW(k\r\)} , (V.77) 



Q-(A- ; r) 4 17-1 



7T 



x + 2 /a; + 2 



(V.78) 



► Quantization of Q modes. 

As in the standard text book material in cosmic perturbation theory, the coefficients 
df may be fixed using the canonical normalization of the modes in the Minkowski, deep 
subhorizon kr — > — oo regime. As discussed, in this limit Q + , which has an oscillatory 
behavior, is the only quantum field (Q_ has an exponentially damping behavior and is 
hence not a quantum mode). We should stress that, of course not all coefficients df and bf 
are fixed by the quantization normalization condition. To fix them, as we will do so below, 
we should impose the constraints ( IV. 5 7ft and (1V.58|) in both superhorizon and asymptotic 



past regimes. Note also that fulfilling these constraints is equivalent to maintaining the 
diffeomorphism and remainder of the gauge symmetry of the system; fluctuations both at 
classical and quantum levels must respect them. 

From the second-order action given in Appendix A, the discussions alluded to above, 
and using the constraints (TV. 591) and (IV. 641) . after some lengthy straightforward algebra the 
momentum conjugate to Q + mode Pq + is obtained as 

Pq + = 2Q' + , (V.79) 

in the kr — > — oo limit, and hence canonically normalized field is \/2Q + . Imposing the usual 
Minkowski vacuum state for v / 2Q+ 

1 



fixes the d~t coefficients 



Q+ ~ = e~ ikl 
2y/k 



4 = -\=, 4 = o. (v.8o) 



The constraint equation flV.61j) can then be used to fix bf coefficients in (TV 741) : 



bf = V2H*l>pt_J± , b+ = 0. (V.81) 
(7 + 2) 



We stress that with the above choice for bf we have only satisfied the constraint f ]V.61j) 



in the asymptotic past kr — > — oo limit. As we will show momentarily, the constraint (IV.61j) 
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in the superhorizon kr — > limit, where Q_, as well as Q+ and assume constant 
values, can only be satisfied when the nonoscillatory Q_ and modes are also taken into 
account. 

The coefficients and b~ , although not fixed by quantization normalization, are related 
through the constraint equation (I V . 6 1 1) . Recalling the exponential growth of the modified 



Bessel function I n (y) for large y [15] . we learn that 

b- = d\ = . 



Demanding the constraint equation flV.6ip to be satisfied in the asymptotic past kr — > — oo 
then leads to 

b 2 ~ J^Lh1><% . (V.82) 

To fix c?2 an d ^2 we should now demand the constraint ( IV. 6 II) to be satisfied in the 
superhorizon kr — > regime too. This will, however, also involve the "oscillatory" Q + and 
ty + modes. Putting these together and making use of the behavior of the Bessel functions 
Z n (y) (where Z n is either H^, or K n ) for small argument y — > [l5], we obtain 

3/4 

y/2 V 3 7 " 



d~2 * "4 ( ^ ) • (V.83) 



► The curvature power spectrum. 

Having fixed all the coefficients bf and df we are now ready to compute the power 
spectrum of metric and curvature perturbations. The power spectrum for the metric per- 
turbations is given by 

which on the large scales (k <C aH) is 



2( 7 + l)(7 + 2) 2 \ n 

As is implicit, we have assumed slow-roll approximation and the value of e and 7 essentially 
remain constant during slow-roll period. The power spectrum of the comoving curvature 
perturbation "Jl, 31 = \& — p ffp , is hence 



2e(7 + l)( 7 + 2) 2 V 7T 
and becomes constant on super-Hubble scales. 
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The spectral index of the curvature perturbations, — 1 = 3 — 2vq, to the leading order 
in the slow-roll parameters is 

37-1 

n»-l~-2^— -e, (V.87) 
7 + 1 

where (1V.71I) and (IV. 761) have been used. We note that the spectral tilt ( IV. 8 711 . depending 
on the value of 7 can be positive or negative; in order to have a red tilt we should consider 
7 > 1/3. 



2. Tensor modes 



To analyze the tensor modes h and T and the action (1V.50|) . it proves useful to decompose 
T into h and a new variable w 

T = Ah + w, (V.88) 

where A is a constant (to be determined). In terms of w the second-order action (1V.5Q|) 
takes the form 

1-2A 



5S. 



(2) 
T 



| / d 3 x (It 



2Ai\) 2 (w'ti + \'Kwh' -k 2 wh + (-f + l)(- 



+ ~(1 + AA 2 ip 2 )(h' 2 - k 2 h 2 + Z -h 2 ) + i/j 2 (w' 2 - k 2 w 2 + "-w 1 )( V.SO) 



where 



— = 'K 2 [2 + 2i\) 2 (1 + 2A - A 2 + (AA - 1)7) - (1 + 8A 2 )e] 

2 



-J{ 2 (2-6)(7 + l) 



From the action (IV. 8911 one can read the equations of motion for h and w 

if / 1 1—94 

w" + (k 2 -^-)w = A( -h" + jXh' -k 2 h+( 1 + 1) (^—f^) X 2 h 

h" + (k 2 - Z ^)h = -yr^r^ (w" + I(5C«,y + ^ - (7 + 1)(^V„ 



(V.90) 
(V.91) 



(V.92) 



(To obtain these results we have used the fact that if) is almost a constant to first order in 
e, cf. ( 1111131) and ( 1III.17I) .) The above equations imply that h and w have both oscillatory 
behavior e lkT in the asymptotic past kr — > —00 region. However, since is negative 

while z"/ 2 is positive, they behave differently in superhorizon kT — > limit; ^ freezes out 
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and — decays. Therefore, in this limit the leading contribution to the right hand side of 
equation of motion for w, which is of order (kr)~ 3 , should vanish. That is, 

(-2 + i)(2 + 7 )^ 2 ^^0, (V.93) 

which implies A=\. This choice for A has an interesting and natural geometric meaning, 
recalling the form of our ansatz for the background gauge field, A a i = ipe a i , where ip is 
a scalar (effective inflaton field) and ef are the 3D triads, and that the triads are "square 
roots" of metric. Perturbing the ansatz and considering only the metric tensor perturbations 



hij, we have 

1 

—i 

2 

Then, recalling (IV. 3 j) and the definition of tij, this implies that A — | naturally removes the 



8e a . = -h^i . (V.94) 



part of the gauge field tensor perturbations which is coming from the perturbation in the 
metric, and hence the "genuine" gauge field tensor perturbation is parameterized by w. 

Inserting T = | + w into the expression (1V.48P and noting the exponential suppression 
of w in the superhorizon scales, 

a 2 rr T ~ , (V.95) 
in the leading order in e and in the kr — > limit. Therefore, the equation of motion for h is 

h" + (k 2 -^)/i~0, ^ = :K 2 (2-e), (V.96) 

which is exactly the same as the standard single scalar field inflation result [3j. The w modes 
do not contribute to the tensor mode power spectrum. 



In the slow-roll approximation flV.96j) is a Bessel equation, whose solution after imposing 
the standard Minkowski vacuum normalization on the kr — > solutions [3], takes the form 



ft = Xlp_ e V + ^)*/* H m {kM) _ (v.97) 

where vt — | + e - Thus, the gravitational waves power spectrum to leading order in slow- roll 
and on super-Hubble scales is set by 

( H\ 2 

y T ~ 2 - \ k=aH . (V.98) 



V 71 

The spectral index of tensor perturbations, tit is given by 



n T ~ -2e . (V.99) 
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To summarize, the gravitational power spectrum and the tensor spectral index of gauge- 
nation are just the 7t and ut of the single scalar model and the tensor-to-scalar ratio r for 
our model is 

gr = 4( 7 + l)( 7 + 2)» 

This result, which can be written as r = — 2 ^^+\)^ nT ^ ma y be contrasted with the usual 
single field consistency relation r = — 8n T [161 ] . 

VI. FITTING GAUGE-FLATION RESULTS WITH THE COSMIC DATA 

We are now ready to confront our model with the observational data. As discussed our 
model for a specific range of its parameters allows for slow-roll inflation and for comparison 
with the observational data we use the results obtained in the slow- roll regime. First, we 
note that in order for inflation to solve the flatness and horizon problems it should have 
lasted for a minimum number of e- folds N e . This amount of course depends on the scale 
of inflation and somewhat to the details of physics after inflation ends |2j. However, for a 
large inflationary scale, like H ~ 10~ 4 — 10 _5 M p i, it is usually demanded that N e ~ 60. As 
a standard benchmark we use N e = 60. 

As for the CMB data, current observations provide values for power spectrum of curvature 
perturbations 7^ and its spectral tilt n% and impose an upper bound on the power spectrum 
of tensor modes 3V, or equivalently an upper bound on tensor-to-scalar ratio r. These values 
vary (mildly) depending on the details of how the data analysis has been carried out. Here we 
use the best estimation of Komatsu et al. jl7] which is based on WMAP 7 years, combined 
with other cosmological data. These values are 

IPs - 2.5 x 10~ 9 , (VI. la) 

n^ = 0.968 ±0.012, (VI. lb) 

r<0.24. (VI.lc) 

Our model has two parameter g and k, and our results for physical observables depend 
also on others parameters which are basically related to the initial values of the fields we 
have in our model. Out of these parameters we choose H, the value of Hubble, and ip, the 
value of the effective inflaton field at the beginning of, or during, slow-roll inflation. The 
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values of other parameters, e, 7 and 5 (initial velocity of the if) field (1III.13j) ). are related to 
these two through flIII.16j) - flIII.19j) . For convenience let us recollect all our results: 



7+1^7+1 (VL2) 
Ze 7 

nx-l--!- (3 7 - 1) In 2+1 , (VI.3) 

T T 1 /2( 7 + l)( 7 + 2)\ 2 l 7 + I 



2N e v 27 + 1 y ' 7 

1 (2 7 +1) 2 # 2 (2 7 + l 



In ^— , (VI.4) 

27r 2 e \M p J (7 + 2) 2 (7 + 1) " 2vr 2 7(7 + 1) 2 ( 7 + ^ ' (VL5) 

Interestingly, the spectral tilt and the number of e-folds only depend on 7 and e and 
hence their values may be fixed using these two, leading to 

0.85 <7< 6.35, e= (0.9- 1.2) x 1(T 2 . (VI.6) 

Notice that lower value for e corresponds to 7 = 6.35. We may now use these values and 
the COBE normalization to read H and i/j: 

jj) = (3.5 -8.0) x KT 2 M pl , H ~ 3.5 x KT 5 M pl . (VI.7) 

(We note that the variation in the value of H over the range (1VI.6I) is about 0.3%.) Finally 
one may use definition of 7 (IIII.19I) to compute g and k 



a 2 

£_= (0.13-5.0) x 10~ 7 , « = (4.6 - 17) x 10 13 MT 4 . (VI.8) 

47T 1 ' 



One can easily check that with the above values for our parameters, our model predicts 
"Pt = 2.45 x 10~ 10 corresponding to r ~ 0.1, which is well within the range to be observed 
by the Planck satellite. 

As mentioned, however, the value of N e is not exactly 60 and may be smaller or larger. 
This possibility has been explored in Fig. ^considering the current observational data (IVI.lb ) 
and (IVI.lb ). The current observational value for — 1 and the bound on r, allows for 
AO < N e < 145 and 0.05 < r < 0.24. In other words, r > 0.05 is one of the predictions of 
our model, which can be tested by the Planck. As handy relations to remember, for 7 = 1, 
r = and for 7 = 10, r = 

It is also instructive to compare our model with usual single field inflationary results. For 
the comparison it is more convenient to replace 7 for 77. Recalling flIII.16j) and flIII.17j) we 
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FIG. 4: The shaded region specifies the range for 7 and iV e allowed by the current observations. 
The dot-dashed line corresponds to (|VI.4|) with r = 0.24. The upper solid curve corresponds to 
(jVL3|) with n s = 0.968 + 0.012 and the lower dashed curve to n K = 0.968 - 0.012. As we see, the 
current observations allows for N e as large as 145. One can show using (|VI.4|) that for the values 
in the shaded region the tensor-to-scalar ratio is in the range 0.05 < r < 0.24. 



learn that 



7 



e — rj 

7] 



(VI.9) 



In our model by definition 7 is a positive number and as is seen from FigJH it is of order 
1 — 10. In terms of e and 77 



1 m. e 



2e(e + r/) 



n 



(VI. 10) 



2r ] '"e-r ] 1 ^ 1 \r](2e-r]) J ' ~ T Tl r/ 2^3 2 1 / 3 ' 
As we see the field ip and its variation during inflation are both small and proportional to rj 
while, as in the standard single field models, N e is inversely proportional to rj (or e). This 
is to be contrasted with the usual results that leads to Lyth bound [l8[ where the inverse 
of slow-roll parameters e or rj are proportional to tp 2 and/or its variation during (slow-roll) 
inflation (Aip) 2 , with ip being the inflaton field. Therefore, our model can naturally produce 
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large tensor-to-scalar ratio, r > 0.01, which will be detectable by Planck satellite, without 
the need for super-Planckian field or field variations. 

VII. SUMMARY, DISCUSSION AND OUTLOOK 

In this work we have presented a detailed analysis of the gauge-flation model which we 



introduced in [l]. We first showed that non-Abelian gauge field theory can provide the 
setting for constructing an isotropic and homogeneous inflationary background. We did so 
by using the global part of the gauge symmetry of the problem and identified the SU(2) 
subgroup of that with the rotation group. We argued that this can be done for any non- 
Abelian gauge group, as any such group has an SU (2) subgroup. Therefore, our discussions 
can open a new venue for building inflationary models, closer to particle physics high energy 
models, where non-Abelian gauge theories have a ubiquitous appearance. 

The Yang-Mills theory cannot serve the job of building inflationary models, and we have 
to consider more complicated gauge theory actions. Among the obvious choices, we have 
checked non-Abelian version of Born-Infeld action 10 (with the symmetric trace prescription 



20(|), which does not lead to a slow- roll dynamics within its space of parameters. We 
have checked F 4 terms which appear in one loop level effective gauge theory action. If we 
parameterize such F A terms as Tr(aF A + f3(F 2 ) 2 ), our analysis shows that it is possible to 
get slow-roll inflationary background for specific range of a and ft parameters. With the 
gauge group SU(2), upon which we have mainly focused in this work, the Tr(F A F) 2 that 
we have considered here can be obtained from specific choices of a and /3. 

As discussed our motivation for considering a Tr(F A F) 2 term was primarily providing 
an explicit, simple realization of our gauge-flation scenario which can lead to a satisfactory 
slow-roll inflation; in this work we were not concerned with explicit derivation or embedding 
of this term from particle physics models. At technical level this happens because the 
dependence of this term on the background metric g^ v appears only through det g and as a 
result the contribution of this term to the energy-momentum of the background will take the 
form of a perfect fluid with P = —p equation of state, perfectly suited for driving an almost 
de Sitter expansion. It is, however, important to study appearance of this K-term through a 



10 



For an analysis of non-Abelian Born-Infled theory within the FRW cosmology see 19 1. 
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rigorous quantum gauge field theory analysis and in particle physics settings. From particle 
physics model building viewpoint, a TriFA F) 2 type term can be argued for, considering 
axions in a non-Abelian gauge theory 2JJ and recalling the axion-gauge field interaction 
term L axion ~ j^TtF AF. Then, integrating out the massive axion field (p leads to an action 
of the form we have considered. If we adopt this point of view our k parameter is then related 



to the cutoff scale A as j| ~ A" 4 |2l|, and hence leading to A ~ lCT 3 Mpi ~ 10 15 GeV. 
In order for this proposal to work, some points should be checked: recalling that H < 10 14 
GeV, A ~ 10H . For this one loop effective action description to make sense it is crucial 
that the cutoff A becomes larger than H, because only axion configurations with subhorizon 
momenta (k > H) will contribute to (quantum) loop corrections. The superhorizon modes, 
as in any quantum field theory on (almost) de Sitter background, are frozen and have become 
classical, and hence do not contribute to quantum corrections. It is also crucial that we are 
in a perturbative regime of the gauge theory with g ~ 10~ 3 . Therefore, we need not worry 
about complications of dealing with a confining (non-Abelian) gauge theory. In our case, 
we are in a weakly coupled regime where the theory is in deconfined phase. We also remark 
that, as argued, during slow-roll inflation regime the contribution of the K-term to the energy 
density of the gauge field configuration should dominate over that of the Yang-Mills part. 
In order for the mechanism for generation of the K-term sketched above to work, one should 
argue how the other possible higher-order terms, at F A level and higher loops (leading to 
higher powers of F in the effective action), are suppressed compared to the K-term. These 
issues will be discussed in a later publication [22]. 

Another interesting feature of our gauge-flation model is its naturalness; that demanding 
to have a successful inflationary model compatible with the current data leads to parameters 
which are within their natural range: the Hubble during inflation H is of order 10 14 GeV, 
and cutoff scale of the theory A ~ 10 15 — 10 16 GeV which are natural within the (SUSY) 
GUT models. Moreover, as is required by the consistency of the theory H is less than cutoff 
A (by one order of magnitude). The other parameter of the theory, the gauge coupling 
g ~ 10~ 3 — 10 -4 , although a bit lower than the value expected for the coupling at the 
gauge unification scale, is also in a natural range. The field value ipi and its displacement 
during inflation ipi — ipf, are both of order 10~ 2 M p \, well within the sub-Planckian regime. 
Therefore, as discussed, the arguments of standard single field inflationary models and the 



Lyth bound 



18] do not apply to our model and we do not face the super-Planckian field 
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problem, which is a generic feature of large-field inflation models, such as chaotic inflation, 
causing concerns about the validity of using classical Einstein gravity. We also note that the 
energy density during inflation 3H 2 M^ ~ (2 x 10 16 GeV) 4 , is the same order as the SUSY 
GUT scale. 

Our other motivation for studying the gauge-flation scenario, which is at least in spirit 
close to beyond standard particle physics model settings, was to provide a setup to address 
cosmological questions after inflation. As we discussed and is also seen from the phase 
diagram in Fig. (TJ after the slow-roll ends we enter a phase where the dynamics of the effective 
inflaton field, and gauge fields in general, is governed by the Yang-Mills term. The effective 
inflaton ip starts an oscillatory phase and through standard (p)reheating arguments, e.g. see 



131 ] . it can lose its energy to the gauge fields. If we have an embedding of our gauge-flation 



scenario into beyond standard models, the energy of these gauge fields will then naturally 
be transferred to all the other standard model particles via standard gauge interactions. 
Therefore, our gauge-flation provides a natural setting for building (p)reheating models, to 
which we hope to return in future works. 
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Appendix A: second-order action 



After a tedious but straightforward calculation, which is also confirmed by the Maple 
codes, the total action to second-order in perturbations is 



.3. Kg 2 (j)\Q 2 . 6. Kg 2 (j)\ g 2 (f) 3 , 4k0 2 Kg 2 (f)\k 2 ^ Q 
- 1 + -if- — + 6^ -1 + -if- * - ifv- 1 + -4- + -if- — aM - 
2 V a 4 J a 2 v a v a 4 ' a 2 a 2 a 4 7 a 2 y a 



fc 2 % 2 2 k0 2 6/^ 2 3 <^Q 2 2^ 2 3 ^ 2 fc 2 • Q 

+ ( — 5- 7T- (1 5-) (^r))-7 n — 5-aJVL- 

a 2 ar cr a cr a z a 2 a a a 



f 12g 2 (j) 3 , «0 2 NT 2g 2 <t> 2 ,Kg 2 (j> 4 3/t0 2 N ^ 2 N Q £ 2 2 _ 
cr a 2 a 2 a 4 cr a 2 a a 1 



o 2 2 o 2 </> 3 l,g 4 (p 5 { k 2 . 2,Kg 2 ^ k0\ £ 2 4 1/^ 2 4 \^k 2 
a 2 a 2 a a0 a 2 3 a 4 a 1 a 2 (p 2 6 a 4 2 a 2 



2# 2 </> 3 / 2/^ 2 ng 2 ^ ^k 2 /( /> 2 o 2 4 ft^W, 

cr cr a 4 a 2 cr a 4 a 4 a 2 



2 



- 3^ 2 + ^ 2 + 6i^ 2 }. (A.l) 
a 2 

Note that the above action is computed fixing the Newtonian gauge E = B = and imposing 
7T, = 0. 



Appendix B: Validity of ansatz (|V.59j) 

In this appendix we prove the validity of ansatz (1V.59|) . In the asymptotic past kr — > —00 
(Minkowski) limit and during the slow-roll inflation, equation (IV.37j) takes the following form 

(Q> - ^(1 + 2 -)g 2 ^ 2 aM) * -f-(Q' + ^(aM)"). (B.l) 

which can be rewritten as 

(0' - ^(1 + ^V«M) ~ ^ 2 Q' - (aM/, (B.2) 

where W is a constant. On the other hand, from (1V".57|) we learn that ?/> should include a 
term ^ = g 2 ip 3 a'M. s , which after combining with (TV. 581) . yields 

(1 + - + -^ 2 )M S = 0, (B.3) 

7 7 

and since 7 > 0, it implies M s = 0, proving the validity of the ansatz (1V.59j) . 
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Appendix C: Short review of constraint structure in cosmic perturbation theory 



Consider the scalar sector of metric perturbations (IV- lft . parameterized as 

ds 2 = -(1 + 2A)dt 2 + a 2 ((1 - 2C]Sij + Id^E) dx i dx j , (C.l) 

where <9j denotes partial derivative respect to x % . For multi-scalar field inflationary models, 
minimally coupled to gravity, generally we have 



Sp = 


i 


Aft) 




5P = 


/ 


Afa) 




6q = 


i 






n s = 


a- 2 (C-A) = 


0, 





(C.2) 



where Vj is equal to jp^ and for simplicity we chose the Newtonian gauge (E = B = 0). 
Also, the gravitational field equations are as below, two constraints 

5q + 2{C + HC) = 0, (C.3) 
k 2 

5p-3H5q + 2—C = 0, (C.4) 
a z 

and a dynamical equation 

5P + 6q + 3H6q + (p + P )C = 0. (C.5) 
From the combination of ( 1C.2j) . (1C.3j) and ( 1C.5|) . we obtain 



£ (0/ + 3H<p! + Vi) = 0, (C.6) 
/ 

which is a summation of the background field equations 

V<pi: <p I + 3H^ I + V I = 0, (C.7) 

That is, (1C.3P is not an independent equation and only the constraint (10.4[) and the dy- 
namical equation (1C.5j) are independent among the Einstein equations. At last in these 
models, the complete set of the equations consists of constraint (]C.4j) . and n field equations 
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of Sifj. This makes it possible to decompose each arbitrary field perturbation into an adia- 
batic (curvature) perturbation along the inflaton trajectory and n — 1 entropy (isocurvature) 
perturbations orthogonal to the inflaton in field space 141 ]. 

In contrast to the multi-scalar field models, gauge-flation has three independent gravita- 
tional field equations: two constraints (1C.3P -( IU~4|) . and a dynamical equation (1C.5P for three 
unknowns (two from the perturbations of the gauge field and one from the fluctuations of the 
metric). In the slow-roll regime, we can write the 2 constrains and the dynamical equation 
in from of (1V.51j) -( TV.52j) and (1V.54|) . So number of the dynamical d.o.f of gauge-flation is 
equal to that of the single scalar field models. Consequently, since in gauge-flation we have a 
non-Abelian gauge field and not simply some scalar fields, the standard approach for dealing 
with multi-scalar fields is not applicable even in the study of scalar sector of gauge-flation 
perturbation theory. 
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